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A  2D  axial-symmetrical  model  is  developed  for  a  button  cell  test  system,  integrating  the  detailed 
electrochemical  elementary  reactions  with  multiple  transport  processes.  The  model  is  validated 
using  experimental  polarization  curves  of  the  cell  under  different  temperature  conditions.  Exten¬ 
sive  simulations  are  performed  to  elucidate  the  complicated  interactions  between  fuel/gas  species, 
adsorbed/desorbed  surface  ions,  elementary  reactions,  and  their  effects  on  the  cell  performance.  Results 
indicate  that  the  concentrations  of  surface  adsorbed  species  of  0(S),  0H(s)  and  H20(S)  ions  presented  at 
the  anode/electrolyte  interface  are  relatively  high,  while  that  of  H(S)  ion  is  relatively  uniform  within  the 
entire  anode.  With  increasing  the  operating  temperature,  the  concentrations  of  surface  adsorbed  0(s), 
0H(S)  and  H20(S)  ions  at  the  anode/electrolyte  interface  are  significantly  improved,  while  the  H(S)  ion  is 
slightly  influenced.  The  adsorbed  surface  species  0(S),  0H(S)  and  H20(S)  ions  are  very  sensitive  to  the  fuel 
compositions  and  cell  voltage. 

Published  by  Elsevier  B.V. 


1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  is  one  of  clean  energy  technolo¬ 
gies  that  can  convert  the  chemical  energy  of  fuels  into  electricity 
in  a  highly  efficient  and  environmentally  benign  manner  [1,2]. 
The  electrochemical  reactions  at  solid/gas  interface  in  porous  elec¬ 
trodes  play  a  very  important  role.  Due  to  very  aggressive  operating 
conditions,  direct  measurement  becomes  difficult  and  the  mod¬ 
eling  technique  has  been  widely  used  to  study  the  complicated 
electrochemical  reaction  processes  [3].  Usually,  a  set  of  conserva¬ 
tion  equations  are  used  to  describe  the  processes  of  heat  transfer, 
fuel/gas  species  transport,  and  charge  (ion  and  electron)  migration. 
These  processes  are  coupled  together  at  solid/gas  interface  through 
Butler-Volmer  equation  [4,5].  In  open  literature,  the  Butler-Volmer 
equation  is  generally  applied  using  bulk  gas  species,  the  detailed 
interim  reaction  species,  which  can  improve  the  underlying  mech¬ 
anism  understanding  of  electrochemical  reactions,  are  usually 
neglected. 

To  overcome  the  limitations  of  the  modeling  method  men¬ 
tioned  above,  elementary  reaction  kinetic  modeling  is  employed 
[6-9].  The  elementary  reaction  models  focus  on  the  description  of 
adsorption/desorption  and  surface  reaction  processes  taking  place 
at  a  local  point,  and  is  able  to  provide  more  precise  electrochem¬ 
ical  charge-transfer  calculation  through  the  balances  of  reaction 
species.  It  is  usually  assumed  that  such  elementary  reactions  are 
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identical  across  the  entire  porous  electrodes.  While  this  assump¬ 
tion  can  simplify  the  modeling  and  related  computations,  it  is  not 
consistent  with  the  fact  that  the  amplitude  of  such  reactions  could 
be  different  from  one  reaction  site  to  another  due  to  the  limitations 
of  mass  transport  processes.  Therefore,  there  is  a  need  to  integrate 
the  elementary  reaction  at  each  individual  site  with  complicated 
fuel/gas  transport  processes  [10-12]. 

The  objective  of  this  paper  is  to  develop  a  mathematical  model 
that  can  link  the  elementary  reactions  at  individual  sites  and  the 
complicated  multiple  transport  processes  in  porous  electrodes.  The 
model  is  validated  using  the  experimental  data  of  a  button  cell. 
Comprehensive  simulation  studies  are  performed  to  investigate 
multi-scale  interactions  of  electrochemical  elementary  reactions 
and  multi-transport  processes. 

2.  Model  setup  and  governing  equations 

Shown  in  Fig.  1  is  the  schematic  illustration  of  a  button  cell  test 
stand.  An  anode-supported  button  cell  Ni-YSZ/YSZ/LSM  is  mounted 
on  one  end  of  a  large  ceramic  tube.  The  humidified  hydrogen  as  a 
fuel  is  supplied  to  the  anode  electrode  via  a  small  ceramic  tube.  The 
surplus  gas  and  product  in  the  anode  flow  out  through  the  large 
ceramic  tube.  The  cathode  electrode  is  exposed  to  the  ambient  air. 
A  test  furnace  is  used  to  control  the  temperature  of  the  cell  test 
stand.  Such  an  experimental  setting  is  employed  as  a  physical  base 
for  mathematical  model  development.  Due  to  the  axial  symmetry, 
two-dimensional  computational  domain  is  considered  to  model 
the  button  cell  test  system.  The  governing  equations  of  the  model 
include  elementary  reactions  and  multiple  transport  processes. 
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Nomenclature 

A  the  Arrhenius  parameter  for  the  rate  constants 

C/  concentration  (mol  m-3 ) 

Dmk  ordinary  diffusion  coefficient  (m2  s-1 ) 

Knudsen  diffusion  coefficient  (m2  s-1 ) 

Dsurf  surface  diffusivity  of  specie  i  (m2  s_1 ) 

Dt  thermal  diffusion  coefficient  (kg  m-1  s_1 ) 

E°  the  standard  potential  (V) 

£act  the  thermal  activation  energy  (kj  mol-1 ) 

F  Faraday’s  constant  (96,485  C  mol-1 ) 

G  the  Gibbs  free  energy  (J) 

ip  Faradic  current  density  (A  m-2) 

kf,kb  the  forward  and  backward  reaction  rate  coefficients 

I<  permeability 

I<e q  the  equivalent  thermal  conductivity  (W  m-1  K-1 ) 

Mf  molecular  weight  of  species  i 

m,n  the  stoichiometric  coefficients 

Pi  the  partial  pressure  of  specie  i 

Qh  the  heat  source  term  (J  mol-1 ) 

r,-  the  rate  of  elementary  reaction  i 

R  universal  gas  constant  (8.314  J  mol  K-1 ) 

Ri  mass  conservation  source  term  of  species  i 

Rn  the  net  reaction  rate 

/tpb  three  phase  boundary  length  per  unit  volume 

(m  m-3) 

T  temperature  (I<) 

V  cell  voltage  (V) 

Xi  molar  fraction  of  specie  i 

z  the  number  of  electrons 

r  the  available  surface  site  density  (mol  cm-2) 
rjact  the  activation  overpotential  (V) 

Oi  surface  coverage  of  adsorbed  specie  i 

p  density  (kg  m-3) 

cre  i  electronic/ionic  charge  conductivity  (S  m-1 ) 

^surf  the  number  of  surface  sites  occupied  by  specie  i 
Wj  mass  fraction  of  species  i 

/jl  viscosity 

£  porosity 

r  tortuosity 

0  volume  fraction  of  electronic  conducting  material 

Subscripts 
An  anode 

Ca  cathode 

El  electrolyte 

Cl  channel 

eff  effective 

Ohm  Ohmic 

elec/chem  Electrochemical/chemical 


2.1.  Elementary  reactions 


The  reactions  in  two  directions  take  place  simultaneously  until 
an  equilibrium  state  is  reached.  The  equilibrium  state  is  dependent 
on  the  local  multiple  processes  and  thus  the  operating  conditions 
of  the  cell. 

At  the  cathode  side,  oxygen  ions  0(Ysz)2-  are  conducted  toward 
the  anode  through  the  oxygen  vacancies  of  electrolyte,  e.g.,  YSZ, 

0(YSZ)X  +  (YSZ)  O(YSZ)2-  +  V(YSZ)#* 

where  0(Ysz)X  is  a  lattice  oxygen  in  YSZ  and  V(YSZ)##  is  an  oxygen 
vacancy  on  YSZ  surface.  At  three  phase  boundary  (TPB)  site,  where 
nickel  and  YSZ  as  well  as  gas  meet  together,  charge  transfer  reac¬ 
tions  take  place  and  several  transfer  pathways  could  occur  at  this 
site.  It  is  recognized  that  incorrect  predictions  may  be  given  by  the 
two-step  charge  transfer  mechanisms  when  the  anode  is  polarized 
at  high  H20  partial  pressures  [14].  Therefore,  in  the  present  ele¬ 
mentary  model  a  single-step  oxygen-spillover  process  is  assumed, 
i.e., 

O(ysz)2  +  Ni  0(s)  +  (YSZ)  +  2e(Nj)- 

Here,  the  oxygen  ions  release  electrons,  forming  oxygen  atoms 
adsorbed  onto  the  surface  of  YSZ.  The  electrons  are  released  to  the 
electronic  conducting  phase  Ni.  When  the  atomic  oxygen  0(s)  on 
the  surface  of  YSZ  meets  the  adsorbed  atomic  hydrogen  H(s)  on  the 
Ni  surface  at  the  TPB  site,  surface  reactions  take  place.  At  this  stage, 
there  exist  three  possible  surface  reactions:  the  atomic  hydrogen 
H(S)  combines  with  atomic  oxygen  0(s),  forming  adsorbed  hydroxyl 
OFI(S) ;  the  hydroxyls  may  combine  together  to  form  water  molecule 
and  atomic  oxygen;  the  OFI(S)  may  further  react  with  atomic  hydro¬ 
gen  H(S)  to  form  adsorbed  water  molecule, 


Electrochemical  reactions  involve  multiple  steps  in  the  anode 
and  strongly  depend  on  the  anode  materials.  For  the  Ni-YSZ 
composite  anode,  it  is  generally  recognized  that  the  following  reac¬ 
tion  processes  are  recognized  when  hydrogen  is  used  as  the  fuel 
[9,13-15]. 

After  hydrogen  diffuses  into  porous  anode,  it  interacts  with 
nickel.  The  hydrogen  can  be  adsorbed  onto  the  surface  of  nickel 
as  hydrogen  atoms;  the  adsorbed  hydrogen  atoms  can  be  desorbed 
off  the  nickel  surface  to  form  hydrogen  molecules,  i.e., 


H(s)  +  0(S)^  OEI(s) 

°H(S)  +  OEI(S)  0(S)  +  EI20(S) 

OH(s)  +  H(s)  O  H20(S) 

In  addition,  the  steam  may  also  experience  the  adsorp¬ 
tion/desorption  processes  on  the  surface  of  Ni, 


H2  +  Ni  ^  2H(s) 


EI20  +  Ni(S)  EI20(S) 


Table  1 

H2/H20  surface  reaction  mechanism  [8,9,12,13]. 
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Ni-surface 

A 

n 

E 

Adsorption/desorptionv  (f,b) 

10-2;  5.593  x  1019 

1 

H2  +  Ni(S)  +  Ni(S)  &  H(S)  +  H(S) 

0;  0 

0:88.12 

2 

Surface  reactions 

H20  +  Ni(S)  ^H20(S) 

1  x  10-1 ;  4.579  x  1012 

0;  0 

0:62.68 

3 

O(s)  +  H(S)  &  OH(S)  +  Ni(S) 

5  x  1022;  2.005  x  1021 

0;  0 

97.9:37.19 

4 

OH(S)  +  H(S)  H20(S)  +  Ni(S) 

3  x  102°;  2.175  x  1021 

0;  0 

42.7:91.36 

5 

OH(S)  +  OH(S)  0(S)  +  H20(s) 

3  x  1021:  5.423  x  1023 

0;  0 

100:209.37 

Ni/YSZ-surface  (charge  transfer  reactions) 

A0 

a 

E 

6 

0(YSZ)X 

+  (YSZ)  ^  O(ysz)2  +  V(ysz) 

1.6  x  1022 

0 

90.9 

7 

O(YSZ)2" 

"  +  Ni(S)  0  0(S)  +  (YSZ)  +  2e(Ni)- 

4.9  x  10“6 

0.5 

c 

Rate  constant  of  Arrhenius  equation  written  as:  k  =  ATn  exp (-E/RT),  the  unit  of  E  is  kj  mol-1,  c  -  calculated  or  estimated  from  references. 


Table  1  summarizes  the  elementary  reactions  in  Ni/YSZ  com¬ 
posite  anode.  The  electrochemical  reactions  in  the  anode  can  be 
represented  with  a  generic  formula  [5], 


Then  the  elementary  kinetic  description  of  charge  transfer  can 
be  formulated  as, 

if,  an  =  HTpb  ZRn  ( 5  ) 

where  z>?an  is  the  source  of  the  faradic  current  derived  from  elec¬ 
trochemical  reactions  in  the  anode;  /Tpb  is  the  volume-specific  TPB 
length  of  the  composite  anode.  Because  the  electroactive  species 
involved  in  the  charge-transfer  reactions  are  surface-adsorbed 
intermediates,  to  study  their  variations  of  spatial  distribution,  the 
surface  coverage  0z-s  are  calculated  through  the  surface  diffusion 
process  [6,12], 

<7surf  CTsur^C- 

-j^Ri  =  -V  •  (D?urfV0,),  0t  =  -ip-i  (6) 

where  or3urf  is  the  number  of  surface  sites  occupied  by  species  i;  r 
is  the  area-specific  density  of  adsorption  site  (mol  cm-2);  D?urf  is 
the  surface  diffusivity  of  species  i;  q  is  the  concentration  of  sur¬ 
face  species  i.  The  elementary  reactions  are  influenced  by  the  local 
species  concentrations  through  surface  diffusion  process,  and  the 
species  concentrations  are  influenced  by  not  only  the  local  ele¬ 
mentary  reactions  but  also  the  complicated  transport  processes  in 
porous  electrodes.  To  study  these  coupling  effects,  the  elementary 
reactions  are  linked  to  macro-scale  transport  processes. 

Since  the  cathode  has  much  simple  elementary  reaction  steps, 
the  equation  from  [9]  is  used  to  describe  the  oxygen  reduction 
process, 


n  m 

^^rpRa*  ^^rrij?Uj  +  ze-  (1 ) 

i  J 

where  Ra  is  the  reactant,  Pu  is  the  product,  rp  and  rrij  are  the  sto¬ 
ichiometric  factors  for  the  reactants  and  products,  respectively,  z 
is  the  number  of  electrons.  Using  Arrhenius  formula,  the  forward 
and  backward  reaction  rate  coefficient  kf  and  kb  can  be  calculated 
as  [7,9], 

kf=ATnexp(^^jexp(a^fj  (2) 

kb  =ATnexp  exp  (-(1  (3) 

where  £act  is  the  thermal  activation  energy,  rj  is  the  electrical  poten¬ 
tial  difference  between  the  reactants  and  products.  The  net  reaction 
rate  Rn  can  be  determined  as, 

Rn  =  fc/JTlRa,]"'  -  *IItp«,P  (4) 

i  j 


if,c  a  —  *ca 


(Po2/Po2)°  25 
1  +(Po2/Po2)°'5 


0.5Friact 

RT 


) 


0.5Friact 

RT 


(7) 


where  a  is  a  free  parameter  used  to  fit  the  model  predictions  with 
experimental  data;  p°2  =  4.9  x  108  exp(-2  x  105/RT);  po2  is  the 
partial  pressure  of  local  oxygen;  pac t  is  the  activation  overpoten¬ 
tial  pact  =  rj-  pe q,  p  =  Ve  -  V{  and  peq  can  be  determined  using  the 
reactants  and  products  of  elementary  reactions, 


AG  RT, 

r]eq~w  +  wln 


fU^\ 
2 _ 

HlRa,]"* 


(8) 


here  AG  is  the  Gibbs  free  energy  of  the  reaction.  Considering  that 
the  hydrogen  electro-oxidation  pathway  is  dominant,  Eq.  (8)  can 
be  simplified  as, 


ro  RT  (  Ph2Pq25 

Peq  =  EU  +  —  In  — - ?- 

2 F  \  Ph2o 


(9) 
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Surface:  Molar  fracbon  of  02 


Max:  0.210 
10.21 

0.21 

0.21 

0.209 

0.209 

0.209 

0.209 

0.209 

0.209 

0.209 


Surface:  Molar  fracbon  of  H20 


Max:  0.232 


Surface:  Molar  fracbon  of  H2  Max:  0.970 


0.96 


10.94 
0.92 
0.9 
0.88 
0.86 

10.84 
0.82 

Min:  0.808 


Fig.  3.  The  molar  fraction  distribution  of  bulk  gas  species  within  electrodes:  (a)  oxygen,  (b)  hydrogen  and  (c)  steam. 


2.2.  Transport  processes 

In  addition  to  electrochemical  elementary  reactions,  SOFCs  also 
include  very  complicated  transport  processes,  such  as  fuel/gas 
transport  in  porous  electrodes,  charge  transport  through  material 
backbone,  as  well  as  heat  transfer.  The  mass  and  charge  trans¬ 
port  processes  are  coupled  together  at  reaction  sites  of  porous 
electrodes  through  elementary  reactions.  To  elucidate  these  com¬ 
plicated  interactions,  the  multi-transport  processes  are  further 
modeled. 


Distance  through  SOFC  axis-symmetric  line(mm) 


(VI 


Fig.  4.  The  distribution  of  local  potential  and  current  along  cell  axis-symmetric  line: 
(a)  local  potential  and  (b)  local  ionic  current. 


2.2.1.  Charge  conservation 

The  charge  conservation  includes  both  ionic  and  electronic 
transport  processes  and  can  be  formulated  using  generic  Ohm’s 
law 

-V  ■  (ufWj)  =  ±tV;  n  =  e,  /  (10) 

where  Vn  is  the  potential;  a^f  is  the  effective  conductivity;  e/i 
represent  electron  and  ion  respectively;  z>  is  the  volume-specific 
faradic  current  in  electrodes.  For  composite  electrodes,  e.g.,  Ni/YSZ, 
the  conductivity  can  be  further  calculated  as  crff  =  0(1  -  £)cre/r , 
aff  =  (1  -  0)(  1  -  e)cri/ r,  here  s  is  the  porosity,  r  is  the  tortuosity, 
and  0  is  the  volume  fraction  of  electronic  conducting  phase  in  the 
electrode. 


2.2.2.  Mass  conservation 

Since  the  multi-species  fuel/gas  are  involved,  the  Stefan- 
Maxwell  equation  is  employed  to  describe  multi-species  transport 
in  both  porous  electrodes  as  well  as  fuel/gas  supply  channel, 


V- 


^Dmffc(VXfc+(Xk_Wk)^r 


+  pv  •  Vwm  =  Rt 


(11) 


where  the  subscripts  m  and  k  stand  for  different  gas  species;  w  the 
mass  fraction;  D^fk  is  the  effective  diffusion  coefficient;  Rt  is  the 
source  term.  The  effective  diffusion  coefficient  D^k  is  determined 
by  combining  molecular  diffusion  in  large  pores  and  Knudsen  dif¬ 
fusion  in  small  pores  comparable  to  the  molecular  mean-free  path, 


D 


eff 

m,k 


(12) 


2.2.3.  Momentum  conservation 

To  describe  the  momentum  conservation  in  the  fuel/gas  supply 
channel  and  porous  electrodes,  the  weakly  compressible  Brinkman 


Table  2 

Reaction  rates  and  source  terms. 
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Reaction  rate 

Reaction  source  term 

Energy  source  term 

ri  =  -  /<i,fa[H(s)]2 

Rh2  =  ~r\ 

^Ca  —  Qohm  +  Qeiec/chem 

r2  =  fc2./[H20]  -  fe2,b[H20CsD] 

Rh2o  —  —r 2 

^E1  =  Qohm 

r3  =  ^3,/[H(s)][0(s)]  -  /<3,b[OH(s)] 

RH(s  )  =  2r:  -  r3  -  r4 

^An  —  Qohm  +  Qeiec/chem 

r4  =  /<4,/[OH(s)][H(s)]  -  /<4,b[H20(s)] 

>3 

O 

JX 

II 

JJi 

1 

1 

KJ 

S3 

n 

II 

O 

r5  =fe5,/[OH(s)]2 -fc5,i,[0(s)][H20(s)] 

^h20(s)  =  ri  -  r4  +  r5 

Table  3 

Boundary  conditions. 

Cathode/channel  interface 

Cathode/electrolyte 

Anode/electrolyte 

Channel/anode 

Fuel  channel 

interface 

interface 

interface 

Ionic  charge  balance 

Insulation 

Continuity 

Continuity 

Insulation 

N/A 

Electronic  charge  balance 

Specified  voltage 

Insulation 

Insulation 

0 

N/A 

Mass 

02/N2  mass  fractions 

Insulation 

Insulation 

Continuity 

H2/H20  mass  fractions 

Momentum 

Pressure 

Wall 

Wall 

Continuity 

Flow  rate,  pressure 

Energy 

Temperature 

Continuity 

Continuity 

Continuity 

Temperature 

equation  and  the  Darcy’s  law  are  used  respectively, 


p(v .  W)v  =  V  • 


-Pi 


+  »7(Vu  +  (Vu)T)  -  -  Kdv)  (V  •  u)I 


+  F 
(13) 


22.4.  Energy  conservation 

The  energy  conservation  in  porous  electrodes  is  described  as, 
-V  •  (KeqVT)  =  Qh  -  CLu  ■  VT  (16) 


--VP  =  v  (14) 

[A 

where  fi  is  the  viscosity  of  fuel/gas,  I<  is  the  permeability  and  can 
be  formulated  as  [16]: 


I<  = 


g3dp 

150(1  -s)2 


(15) 


Here  dp  is  the  average  particle  diameter  of  porous  electrodes. 


where  I<eq  =  fiKL  +  (1  -  is  the  equivalent  thermal  conductiv¬ 
ity,  I<l  is  the  effective  heat  conductivity  of  liquid  phase,  KP  is  the  heat 
conductivity  of  solid  phase,  j3  is  the  volume  fraction  of  liquid  phase; 
Ci  is  the  effective  volumetric  heat  capacity;  Qh  is  the  general  heat 
source  term  and  can  be  described  as  Qh  =  Qeiec/chem  +  Qohrm  which 
includes  the  chemical/electrochemical  reaction  heat  (Qeiec/chem) 
and  ohmic  resistance  induced  heat  generation  (Qohm  =  -i  rj  = 
i2R).  Table  2  summarizes  the  relevant  heat  generation  mechanisms 
and  the  corresponding  mathematical  expressions  according  to  the 
equilibrium  of  chemical/electrochemical  reactions. 


Max:  5.919e-9 


:  5.806 
5.692 
5.579 
5.465 
5.352 
5.238 
_5.125 
“5.011 
J4.898 
Min:  4.841e-9 


Fig.  5.  The  concentration  distribution  of  adsorbed  surface  species  within  the  anode:  (a)  0(S),  (b)  H(s)l  (c)  OH(s),  and  (d)  H20(S). 
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Distance  through  SOFC  axis-symmetric  line  (mm) 


Distance  through  SOFC  axis-symmetric  line  (mm) 

Fig.  6.  Effects  of  operating  temperature  on  the  bulk  fuel  species  distribution  along  cell  axis-symmetric  line:  (a)  hydrogen,  (b)  steam,  and  (c)  oxygen. 


2.3.  Boundary  conditions 

To  solve  the  equations  of  elementary  reactions  and  multi¬ 
physics  transport  processes,  boundary  conditions  are  needed. 
According  to  the  experimental  setup,  the  boundary  conditions  are 
summarized  in  Table  3. 

3.  Solution  algorithm  and  model  validation 

The  mathematical  model  is  solved  using  COMSOL  MULTIPH- 
SICS  V4.0.  The  model  parameters  are  provided  in  Table  4.  For  a 
specified  cell  voltage  at  the  cathode  electrode  boundary,  the  dis¬ 
tributions  of  various  parameters  are  calculated.  Accordingly,  the 
average  cell  current  density  is  obtained.  The  cell  polarization  curve 
is  then  obtained  by  specifying  a  series  of  cell  voltages  and  calculat¬ 
ing  the  corresponding  average  cell  current  density.  And  the  species 
associated  with  multi-physicochemical  processes  are  calculated  at 
each  of  the  specified  voltage  conditions. 

The  model  is  validated  using  experimental  polarization  curves, 
where  the  humidified  hydrogen  (97%  hydrogen  and  3%  vapor)  is 


Distance  through  SOFC  axis-symmetric  line  (mm) 


used  as  the  fuel  and  ambient  oxygen  is  used  as  the  oxidant.  The  cell 
polarization  curve  is  obtained  at  temperature  600  °C,  650  °C  and 
700  °C  respectively.  As  shown  in  Fig.  2,  the  model  predictions  match 
with  experimental  data  reasonably  well.  The  validated  model  is 
then  used  for  further  simulation  studies. 

4.  Results  and  discussion 

4. 1.  Distributions  of  fuel/gas  and  local  potential! current 

Fig.  3  shows  the  molar  fraction  distributions  of  hydrogen, 
vapor,  and  oxygen,  where  the  furnace  temperature  is  973  K  and 
the  cell  voltage  is  0.5  V.  The  molar  fraction  of  hydrogen  shows 
little  variation  in  fuel  supply  tube,  but  decreases  in  the  anode 
toward  the  anode/electrolyte  interface.  The  variation  is  well  aligned 
with  the  size  of  the  cathode.  The  molar  fraction  of  vapor  shows 
opposite  trend  and  increases  toward  the  anode/electrolyte  inter¬ 
face.  The  resulting  such  distributions  are  due  to  the  fact  that 
the  electrochemical  reactions  in  the  anode  consume  hydrogen 
while  producing  water  molecules.  Similarly,  on  the  cathode  side, 


Distrance  through  SOFC  axis-symmetric  line  (mm) 


Fig.  7.  Effects  of  operating  temperature  on  the  adsorbed  species  along  cell  axis-symmetric  line:  (a)  0(S),  (b)  H(s),  (c)  OH(s),  and  (d)  H20(S). 
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Fig.  8.  Effects  of  fuel  compositions  on  the  cell  performance. 


the  molar  fraction  of  oxygen  shows  a  slight  decrease  toward  the 
cathode/electrolyte  interface,  in  which  the  oxygen  is  consumed  by 
electrochemical  reactions. 

The  electrical/ionic  potentials  and  ionic  current  distribution  are 
examined  along  the  axis-symmetric  line.  As  shown  in  Fig.  4a,  when 
the  cell  voltage  is  set  at  0.5  V,  the  electrical  potential  of  the  cathode 
and  anode  remains  at  0.5  V  and  0  V  respectively  and  does  not  show 
any  gradient  along  the  axial-symmetric  line.  The  ionic  potential  is 
about  0.5  V  in  the  anode  and  approaches  to  zero  in  the  cathode. 
Obvious  gradients  of  ionic  potential  exist  across  the  electrolyte. 
The  current  distribution  along  the  axial-symmetric  line  is  shown  in 
Fig.  4b.  Clearly  the  ionic  current  reaches  the  highest  value  within 
the  electrolyte,  which  is  consistent  with  ionic  potential  distribu¬ 
tions.  This  observation  suggests  that  the  electrochemical  reactions 
mainly  occur  near  the  electrode/electrolyte  interface.  When  the 
specified  cell  voltage  increases  to  1 V,  the  ionic  current  significantly 
decreases  as  shown  in  Fig.  4b. 


Distance  through  SOFC  axis-symmetric  line  (mm) 


Table  4 

Operation  parameters  [8,17,18]. 


Value 

Units 

Cathode  layer  thickness,  radius 

4  x  10-5,  2.5  x  10~3 

m 

Electrolyte  layer  thickness,  radius 

2  x  10-5,  6.5  x  1(T3 

m 

Anode  layer  thickness,  radius 

3.5  x  10“4,  6.5  x  1(T3 

m 

Inlet  fuel  channel  radius,  wall 
radius 

1  x  10“3, 1.5  x  1(T3 

m 

Outlet  fuel  channel  radius 

6.5  x  10“3 

m 

Faraday  constant 

96485 

C  mol”1 

Gas  constant 

8.314 

J  mol-1  K"1 

Operating  temperature 

700 

°C 

Pressure 

1 

atm 

Anode  ionic  conductivity 

3.34 x 104 exp(-10300/T) 

S  m1 

Cathode  ionic  conductivity 

3.34 x 104 exp(-10300/T) 

S  m1 

Anode  electronic  conductivity 

2  x  106 

S  itT1 

Cathode  electronic  conductivity 

42  x 106 exp(-1150/T)/T 

S  itT1 

Electrolyte  ionic  conductivity 

3.34 x  104  exp(-10300/T) 

S  itT1 

Porosity  (e) 

0.4 

1 

Particle  diameter 

2  x  10“6 

m 

Tortuosity 

[(3  -  £)/2]0'5 

1 

Permeability 

s3d2/ 150(1  -e)2 

1 

Inlet  velocity 

0.2 

m  s-1 

Surface  site  density  of  Ni 

6.1  x  10“9 

mol  cm-2 

Specific  three-phase  boundary 
length 

oo 

X 

m  m3 

4.2.  Surface  species  distribution  within  the  anode 

The  adsorbed  surface  species  such  as  0(s),  H(S),  OH(s),  FI20(S) 
are  very  important  intermediates  influencing  the  electrochemical 
reactions  at  TPB  sites  in  the  anode.  As  shown  in  Fig.  5,  the  concentra¬ 
tion  of  species  H(S)  increases  from  the  electrolyte/anode  interface 
toward  the  anode/channel  interface,  while  the  concentrations  of 
other  species,  e.g.,  0(s),  OFI(S),  H20(s),  decrease.  One  may  also  notice 
that  the  concentration  of  species  H(S)  is  relatively  uniform  within 
the  anode,  however,  those  of  the  species  0(s),  OFI(S),  FI20(s)  show 
relatively  large  gradients  with  the  highest  concentrations  at  the 
electrolyte/anode  interface.  It  is  generally  recognized  that  Ni  is  a 
very  active  catalyst  and  uniformly  distributed  within  the  composite 
anode.  As  a  result,  the  H(S)  ions  are  generated  throughout  the  entire 
anode,  leading  to  the  relatively  uniform  distribution  of  species  H(S). 
The  species  0(s)  is  generated  from  0(Ysz)2-  i°ns  conducted  from  the 
cathode  side.  The  species  0(s)  is  consumed  at  the  anode/electrolyte 
interface  and  keeps  being  consumed  within  the  anode,  resulting 
in  the  decrease  of  spatial  concentration  from  the  anode/electrolyte 


Distrance  through  SOFC  axis-symmetric  line  (mm) 


Fig.  9.  Effects  of  fuel  composition  on  the  adsorbed  species  along  cell  axis-symmetric  line:  (a)  0(S),  (b)  H(s),  (c)  OH(s),  and  (d)  H20(s). 
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Fig.  10.  Sensitivity  of  species  and  cell  parameters  with  respect  to  fuel  compo¬ 
sitions.  The  parameter  values  are  recorded  at  the  intersecting  point  between 
electrode/electrolyte  interface  and  cell  axial-symmetric  line. 


corresponding  species  transport  processes  under  different  temper¬ 
atures,  the  modeling  analysis  is  further  performed.  Fig.  6  shows 
the  distributions  of  fuel/gas  along  the  axial-symmetrical  line  of 
the  cell  at  temperature  873  K,  923  K,  973  K,  respectively.  Essen¬ 
tially,  the  molar  fractions  of  hydrogen  and  oxygen  decrease  from 
the  channel/electrode  interface  toward  the  electrolyte/electrode 
interface.  The  gradients  of  hydrogen  and  oxygen  molar  fractions 
increase  with  increasing  operating  temperature,  indicating  that 
more  hydrogen  and  oxygen  are  consumed  in  higher  operating  tem¬ 
peratures.  On  the  other  hand,  the  concentration  of  vapor  (H20) 
shows  opposite  trend,  and  the  production  of  H20  increases  with 
increasing  the  operating  temperature. 

Fig.  7  shows  the  distributions  of  adsorbed  surface  species  within 
the  anode  along  the  axial-symmetrical  line  of  the  cell  at  temper¬ 
ature  873  K,  923  K,  973  K,  respectively.  The  species  H(s)  ion  shows 
relatively  uniform  distribution  and  a  slight  variation  under  differ¬ 
ent  temperatures.  The  concentrations  of  other  three  species  0(s) 
ion,  OH(S)  ion,  and  H20(s)  ion  are  very  low  but  increase  significantly 
approaching  the  anode/electrolyte  interface,  and  increase  with 
increasing  the  operating  temperatures.  These  observations  again 
indicate  that  the  electrochemical  reactions  in  the  anode  mainly  take 
place  near  the  anode/electrolyte  interface,  even  though  the  com¬ 
posite  anode  is  employed.  One  also  can  see  that  the  concentrations 
of  involved  OH(s)  ions  and  H20(S)  ions  are  much  lower  than  those 
of  H(S)  ions  and  0(s)  ions. 


interface  toward  the  anode/channel  interface.  The  OH(s)  ions  and 
H20(S)  ions  have  much  lower  concentrations  compared  with  those 
of  0(S)  and  H(s)  ions,  and  their  highest  concentrations  occur  in  the 
area  close  to  the  anode/electrolyte  interface.  These  observations 
suggest  that  most  of  the  surface  reactions  take  place  around  the 
anode/electrolyte  interface. 

4.3.  Effects  of  operating  temperatures 

Experimental  results  show  that  the  cell  performance  is  signifi¬ 
cantly  affected  by  the  operating  temperatures  (Fig.  3).  To  study  the 


4.4.  Effects  of  fuel  compositions 

The  effects  of  fuel  composition  on  SOFC  performance  have  been 
widely  investigated  through  modeling  techniques  in  open  litera¬ 
ture,  however  the  elementary  reactions  are  generally  neglected.  As 
a  result,  the  adsorbed  surface  intermediates  are  not  considered.  In 
this  section,  the  surface  electrochemistry  details  are  further  stud¬ 
ied.  Fig.  8  shows  the  polarization  performance  of  the  cell  under 
different  fuel  composition  conditions.  Obviously,  larger  hydrogen 
fraction  in  fuel  supply  leads  to  better  cell  performance.  Fig.  9 


CV=0.2V  CV=0.5V  CV=0.8V 


-0.05  0  0.05  -0.05  0  0.05  -0.05  0  0.05 

Relative  difference  of  parameters  to  +0.01  V  change  on  cell  voltage(CV) 


Fig.  11.  Sensitivity  of  species  and  cell  parameters  with  respect  to  cell  voltages.  The  parameter  values  are  recorded  at  the  intersecting  point  between  electrode/electrolyte 
interface  and  cell  axial-symmetric  line. 


Y.  Xie,  X.  Xue  /  Journal  of  Power  Sources  209  (2012)  81-89 


89 


shows  the  corresponding  distributions  of  adsorbates  along  the 
axial-symmetrical  line  of  the  cell.  Basically,  the  concentration  of 
adsorbed  H(S)  ion  increases  with  increasing  hydrogen  fraction  in 
the  fuel,  and  are  relatively  uniform.  The  concentration  of  adsorbed 
species  0(s)  increases  near  the  anode/electrolyte  interface  with 
increasing  hydrogen  fraction  in  the  fuel,  while  those  of  adsorbed 
species  0H(S)  ion  and  H20(S)  decrease.  When  hydrogen  fraction  in 
the  fuel  is  low,  less  H(S)  ions  are  generated.  This  leads  to  low  chem¬ 
ical  potential  for  0(s)  production.  As  a  result,  the  electrochemical 
reaction  rate  reduces  and  the  cell  performance  decreases.  On  the 
other  hand,  when  more  H20(S)  and  OH(s)  ions  are  generated  by 
higher  bulk  H20  concentration,  it  can  raise  the  backward  rate  of 
surface  chemical  reactions.  Accordingly  the  electrochemical  reac¬ 
tion  rate  will  also  reduce.  Therefore,  the  effect  of  fuel  composition 
on  the  cell  performance  depends  on  not  only  the  reactant  fuels  but 
also  the  product  species  generated. 

To  study  the  sensitivity  of  the  generated  intermediate  species 
and  other  cell  parameters  with  respect  to  the  variations  of  sup¬ 
plied  hydrogen  and  oxygen,  the  simulations  are  run  by  adding  1% 
more  hydrogen  and  oxygen  to  the  anode  and  cathode  respectively. 
Fig.  10  shows  the  corresponding  variations  of  various  parameters 
at  the  anode/electrolyte  interface.  It  is  clearly  seen  that  the  surface 
adsorbed  species  0(s),  H20(s),  and  OH(s)  ions  are  very  sensitive  to  the 
composition  changes  of  fuel/gas.  This  result  indicates  that  increas¬ 
ing  fuel/gas  concentrations  can  effectively  improve  the  chemical 
potentials  for  the  surface  intermediate  species  productions  and 
surface  electrochemical  reactions. 

45.  Effects  of  cell  voltage 

The  sensitivities  of  various  parameters  with  respect  to  the  vari¬ 
ation  of  cell  voltage  are  also  studied.  Fig.  11  shows  parameter 
variations  at  the  anode/electrolyte  interface  when  the  cell  voltage 
increases  from  0.2  V,  0.5  V,  and  0.8  V  to  0.21  V,  0.51  V,  and  0.81  V 
respectively.  The  cell  current  and  three  adsorbed  surface  species 
0(S),  FI20(S),  and  OFI(S)  ions  demonstrate  relatively  high  sensitivi¬ 
ties.  The  concentrations  of  adsorbed  surface  species  H(s)  and  the 
bulk  H2  and  H20  show  slight  variations.  The  variations  of  bulk  02 
and  the  cell  temperature  are  negligible.  One  may  also  see  that  rel¬ 
atively  high  sensitivity  is  obtained  at  cell  voltage  0.8  V,  particularly 
for  adsorbed  surface  species  OFI(S)  and  H20(S).  At  high  cell  volt¬ 
age  conditions,  the  consumed  hydrogen  and  oxygen  are  low.  As  a 
result,  their  concentrations  are  maintained  at  relatively  high  level. 
Accordingly,  high  chemical  potentials  are  generated  for  adsorbate 
productions. 

5.  Conclusion 

A  2D-axial-symmetrical  SOFC  model  is  developed  using  an 
anode-supported  Ni-YSZ/YSZ/LSM  button  cell  test  system  as  the 


physical  base,  in  which  the  humidified  hydrogen  is  used  as  the  fuel 
and  ambient  air  as  the  oxidant.  The  model  considers  the  elementary 
reactions  at  TPB  sites  in  the  anode  and  multiple  transport  processes. 
The  elementary  reactions  are  linked  with  multi-transport  processes 
in  a  coherent  way.  The  model  is  validated  with  experimental  data. 
Simulation  results  indicate  that  very  complicated  interactions  exist 
among  fuel/gas  transport  in  porous  electrodes,  adsorbed/desorbed 
species  on  the  surface  of  porous  electrodes,  and  charge  transport 
through  solid  backbone,  as  well  as  elementary  reactions.  The  con¬ 
centrations  of  surface  adsorbed  species  of  0(s),  OFI(S)  and  H20(S) 
ions  presented  at  the  anode/electrolyte  interface  are  relatively  high, 
while  that  of  H(S)  ion  is  relatively  uniform  within  the  entire  anode. 
With  increasing  the  operating  temperature,  the  concentrations  of 
surface  adsorbed  0(s),  OFI(S)  and  H20(S)  ions  at  the  anode/electrolyte 
interface  are  significantly  improved,  while  that  of  FI(S)  ion  is  slightly 
influenced.  The  adsorbed  surface  species  0(s),  OFI(S)  and  H20(S)  ions 
are  very  sensitive  to  the  variations  of  the  supplied  hydrogen  and 
oxygen  as  well  as  the  cell  voltage. 
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